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Abstract Wc consider the physical origin of the hemispheric pattern of filament 
chirality on the Sun. Our 3D simulations of the coronal field evolution over a 
period of 6 months, based on photospheric magnetic measurements, were previ- 
ously shown to be highly successful at reproducing observed filament chiralities. 
In this paper we identify and describe the physical mechanisms responsible for 
this success. The key mechanisms arc found to be (1) differential rotation of 
north-south polarity inversion lines, (2) the shape of bipolar active regions, 
and (3) evolution of skew over a period of many days. As on the real Sun, the 
hemispheric pattern in our simulations holds in a statistical sense. Exceptions 
arise naturally for filaments in certain locations relative to bipolar active regions, 
or from interactions between a number of active regions. 
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nctic 

1. Introduction 

Solar filaments (also known as prominences) are classified as having either dex- 
tral or sinistral chirality (handedness). The chirality depends on whether the 
axial magnetic field in the filament points right or left, when seen from the 
side of the filament with positive magnetic polarity in the photosphere. Ob- 
servations show that quiescent filaments in the northern hemisphere tend to 
be dextral, while those in the southern hemisphere tend to be sinistral (Rust, 
1967; Leroy, Bommier, and Sahal-Brechot, 1983; Martin, Bilimoria, and Tra- 
cadas, 1994; Pevtsov, Balasubramaniam, and Rogers, 2003). The latter paper 
found that 80% - 85% of quiescent filaments follow this hemispheric pattern, 
which remains the same when the Sun's polar field reverses every 11 years. 

Theories for the formation of sheared magnetic fields in filaments are di- 
vided between those invoking a sub-surface origin and those requiring purely 
surface effects. Rust and Kumar (1994) propose that the helical magnetic fields 
in filaments emerge from the convection zone as already twisted structures, 
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as in the simulations of Gibson et al. (2004). Indeed, observations show that 
emerging magnetic flux may appear at the solar surface carrying electric current 
(Leka et al., 1996), and in a twisted flux-rope geometry (Lites, 2005). However, 
filaments observed so far in such circumstances have been short-lived active 
region filaments, rather than large stable quiescent filaments. For the latter, the 
emerging flux rope hypothesis is inconsistent with the tendency for filaments to 
form between multiple active regions, rather than in centres of flux emergence 
activity (Mackay et al., 2008). 

The alternative origin of sheared fields in filament channels is surface motions; 
we know that after active regions emerge they are advected by the large-scale sur- 
face motions of differential rotation and meridional flow, and the flux is dispersed 
by the action of supergranular convection. The basic principle for formation 
of axial magnetic fields along polarity inversion lines (PILs) was proposed by 
van Ballegooijen and Martens (1989), based on flux cancellation at the PIL 
following surface shearing. More recent models have considered the interaction 
of two bipolar active regions, where the filament forms above a PIL between 
them (Martens and Zwaan, 2001; Litvinenko and Wheatland, 2005; Welsch, 
Devore, and Antiochos, 2005; Mackay and van Ballegooijen, 2005). The latter 
paper allowed the bipoles to have a non-zero magnetic helicity (twist) , effectively 
including the emergence of twisted fields from below the surface. As the helicity 
of active regions is observed to follow a hemispheric pattern (Pevtsov, Canfield, 
and Metcalf, 1995), this hemispheric pattern of twist was suggested to be a 
contributory factor to the hemispheric pattern in filament chirality. 

The work of Mackay and van Ballegooijen (2005), who considered the basic 
interaction between a pair of magnetic bipoles, is the motivation for the present 
study. We have constructed a new version of the numerical model for the evo- 
lution of the global solar corona over multiple Carrington rotations, based on 
observations of photospheric magnetic field distributions. The model has been 
described in detail in Yeates, Mackay, and van Ballegooijen (2007, hereafter 
"Paper 1") and Yeates, Mackay, and van Ballegooijen (2008a, hereafter "Paper 
2"), and is summarised very briefly in Section 2 below. In the second paper, 
we compared the chiralities of 109 observed filaments with the skew angle of 
our simulated coronal field at the observed locations of the filaments. We found 
that the model was highly successful in reproducing the chirality of observed 
filaments, including exceptions to the hemispheric pattern. In this paper we 
consider the specific mechanisms responsible for the chirality in our simulation. 
We have identified eight key mechanisms which produce skewed coronal fields 
above PILs, and each mechanism is described in Section 3 before considering 
their relative importance in Section 4. We conclude in Section 5 with a summary 
of the cause of the hemispheric pattern in our simulation. 

2. Overview of the Model 

Our model for the 3D coronal field evolution is based on that of van Ballegooijen, 
Priest, and Mackay (2000), and has been described in Papers 1 and 2. The mag- 
netic field in the corona evolves over time in response to (1) the emergence of new 
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Figure 1. Surface radial field on day 229, and chirality of simulated filaments during Car- 
rington Rotation 1952. Shading shows the sign of radial magnetic field on the photosphere 
(white for positive, grey for negative), and thin lines show contours of the same field (solid for 
positive, dashed for negative). Coloured lines show simulated skew at the locations of observed 
filaments in the comparison (red for dextral, blue for sinistral, and green for weak skew). 

magnetic bipoles from below the photosphere, and (2) advection by large-scale 
motions on the lower (photospheric) boundary. The newly-emerging bipoles take 
a simple mathematical form and, in order that the simulation remains accurate 
to the observed fields on the Sun, the bipole properties (latitude, longitude, 
size, and magnetic flux) are determined from observations. Figure 1 shows a 
snapshot of the photospheric field on day of year 229, mid-way through the 
6-month simulation. Visible are newly-emerged bipoles at low latitudes, and 
also older regions which have been sheared by differential rotation, have been 
transported poleward by the meridional circulation, and have interacted with 
the surrounding field by surface diffusion (representing the cumulative effect of 
supcrgranular convection) and reconnection. 

The coronal field associated with the photospheric field is evolved using the 
magnetofrictional relaxation technique, effectively through a continuous sequence| 
of nonlinear forcc-frcc equilibria that allows the build up of sheared fields and 
magnetic hclicity. The hclicity in our model originates from three distinct sources 
(illustrated in Yeates, Mackay, and van Ballegooijen, 2008b): (1) new bipoles 
emerge twisted, (2) currents develop at the interface between old and new flux 
regions, and (3) surface motions shear the coronal magnetic field. Note that the 
subsurface origin of emerging hclicity is outside the scope of our simulations, and 
is still uncertain. Theories include the interaction of rising flux tubes with the 
large-scale poloidal field (Choudhuri, Chatterjee, and Nandy, 2004), with helical 
turbulence (the S-effect, Longcope, Fisher, and Pevtsov, 1998), or with the 
Coriolis force (Holder et at, 2004). Observations of current helicity (e.g. Pevtsov, 
Canficld, and Metcalf, 1995; Hagino and Sakurai, 2004) and of X-ray sigmoids 
(Canfield and Pevtsov, 1999) show that the majority twist of active regions is 
negative in the northern hemisphere and positive in the southern hemisphere. 
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Figure 2. Mechanisms which always produce the same type of skew (i.e., dextral or sinistral) 
in each hemisphere. In each example, shading shows the sign of radial field B r on the photo- 
sphere (white for positive, grey for negative), thin lines show contours of the same field (solid 
for positive, dashed for negative), and thick blue lines show selected coronal field lines. 

As described in Paper 2, we have run simulations where the sign of cmerging- 
bipole twist (/?) matches this majority observed sign in each hemisphere, and 
simulations where it takes the opposite sign. The filament hemispheric pattern 
is best reproduced for the majority observed sign of (3. 

The coloured lines in Figure 1 show the skew type (i.e., dextral, sinistral, or 
weak) of our simulated coronal field at the locations which were compared with 
observed filaments during a particular 27-day period. This shows results from 
the simulation run with correct hemispheric sign of bipole twist. The majority of 
these filament locations match the observed chirality, including two exceptions 
to the hemispheric pattern, labelled El and E2 on the figure. The only filaments 
where our simulation produces the wrong chirality are one on the polar crown 
(XI), and another at lower latitude (X2) which is only weakly skewed. In this pa- 
per we address the question of why our simulation is so successful at reproducing 
the observed chiralitics. 

3. Physical Mechanisms Producing Skew 

Among our 109 simulated filaments we have identified eight different mecha- 
nisms which produce skewed coronal magnetic fields above PILs. Three of these 
mechanisms always produce the same type of skew (i.e., dextral or sinistral) 
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Figure 3. Mechanisms which can produce either dextral or sinistral skew depending on in- 
dividual circumstances. In each example, shading shows the sign of radial field B r on the 
photosphere (white for positive, grey for negative), thin lines show contours of the same field 
(solid for positive, dashed for negative), and thick blue lines show selected coronal field lines. 



in each hemisphere, and are illustrated in Figure 2. The other five mechanisms 
may produce cither type of skew depending on individual circumstances, and are 
illustrated in Figure 3. Note that changing the polarity of the magnetic field in 
each figure does not affect the skew type. We describe each mechanism in more 
detail below. 
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Figure 4. Schematic showing how surface diffusion converges field line foot-points towards a 
PIL, leading to an increase in the skew angle. 

3.1. Differential Rotation 

The most conceptually simple way to produce a skewed coronal field is through 
differential rotation of a pre-existing field. The type of skew produced depends 
on the orientation of the initial field and the PIL. Consider first a coronal field 
oriented in the east-west direction, over a north-south PIL. This is mechanism 1, 
illustrated by a cartoon in Figure 2 (for the northern hemisphere) . Here differen- 
tial rotation has little effect on the higher coronal field lines, as both footpoints 
are at the same latitude. However, the photosphcric PIL itself is sheared and so 
the coronal field develops a skew relative to it. Because the equator rotates faster 
than the poles, this skew will always be dextral in the northern hemisphere and 
sinistral in the southern hemisphere, in accordance with the hemispheric pattern. 
This effect of differential rotation of a north-south PIL was noted by Zirker et 
al. (1997), but earlier authors had dismissed differential rotation as a potential 
source of the hemispheric pattern because they had considered a north-south 
field across an east-west PIL (e.g., van Ballegooijen and Martens, 1990; Rust 
and Kumar, 1994; Priest, van Ballegooijen, and Mackay, 1996). This alternative 
orientation of field and PIL is mechanism 2 in Figure 2. In this case the PIL 
is unaffected by the differential rotation, but the footpoints of the coronal field 
lines arc shifted, resulting in the wrong type of skew, sinistral in the northern 
hemisphere and dextral in the southern hemisphere. Thus mechanism 2 always 
acts to produce exceptions to the hemispheric pattern. 

We remark here that the skew developed by differential rotation in mech- 
anisms 1 and 2 (or indeed by any of the other mechanisms described in this 
paper) is enhanced over time by convergence of flux towards the PIL, owing 
to supergranular diffusion. This convergence, sketched in Figure 4, moves the 
footpoints of coronal field lines towards the PIL irrespective of its north-south or 
east-west orientation, thus enhancing any existing skew, although not producing 
it from scratch. The process is well illustrated by the Ha sequence in Figure 4 of 
Wang and Muglach (2007). Note that our simulations do not include any extra 
converging flow towards the PIL over and above that due to diffusion; although 
such localised flows have been invoked in some models of filament formation 
[e.g., Kuijpers, 1997) they are not supported by observations (Hindman, Haber, 
and Toomre, 2006). 
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3.2. Filaments on the Internal PIL of a Bipolc 

In our simulation there are 32 filaments (out of 109) which form internal to a 
single active region. The skew of these filaments shows a characteristic behaviour 
(mechanism 3, Figure 2). Initially the skew depends only on the sign of j3 - the 
bipole twist for the region, but after 1 to 2 Carrington rotations the effects of 
differential rotation lead to dextral skew in the northern hemisphere and sinistral 
in the southern hemisphere, irrespective of the sign of (3. This is essentially 
mechanism 1 in action, but it is assisted by the tilt angle of the bipolc, which 
enhances the latitudinal shearing. So for the majority hemispheric sign of (3 
such internal filaments will always follow the hemispheric pattern, and even for 
the minority sign of (3 the hemispheric pattern will be recovered after 1 to 2 
Carrington rotations (30-60 days). 

3.3. Filaments on the External PIL of a Bipolc 

Many filaments occur along the outside edge (periphery) of a single bipolar 
region, and in the simulation their skew is directly affected by the shape of field 
lines in the bipolc (mechanism 4, Figure 3). This idea is supported by the recent 
observations of filament formation by Wang and Muglach (2007). They show an 
example where Ha fibrils on the outside edge of a newly-emerged active region 
become rapidly aligned with the horizontal field component of the bipole, which 
lies primarily along the bipole's external PIL (see their Figure 5). 

The skew of such a filament depends primarily on where along the bipolc 
boundary the filament lies, and on the polarity of the background field relative 
to the bipole. The field lines are most strongly skewed at the north and south 
ends of the bipole, with weaker skew on the middle part of the external PIL (at 
least until twist has diffused out from the bipole centre). This is seen in Figure 
5, which shows the shape of the bipole for the six possible sign combinations of 
background field and bipole twist (3. The type of skew (dextral, sinistral, or weak) 
is indicated by the colouring along the PIL in each case. Notice that changing 
the sign of (3 for the bipole does not change the type of skew on the external 
PIL, but merely the strength of skew. In our simulation there arc 31 filaments 
occuring at either the north or south end of such an external PIL, and the skew 
pattern in Figure 5 is seen to hold for the majority (24) of such filaments. 

How docs this relate to the hemispheric pattern? In the northern hemisphere, 
the observed hemispheric pattern of hclicity gives (3 < 0. If the background field 
has the following polarity then Figure 5(c) shows that a filament occuring at the 
north end of the external PIL will have weak or slightly sinistral skew, whereas 
one occuring at the south end will have dextral skew. If the background field 
has the leading polarity, then Figure 5(f) shows that a filament at the north 
end will be dextral, whereas at the south end the field is weakly or slightly 
sinistrally skewed. Overall, then, there will be more dextral skew produced 
by this mechanism than sinistral skew in the northern hemisphere. A similar 
argument shows that in the southern hemisphere, where (3 > 0, there will be 
more sinistral skew produced. Thus filaments on the external PIL of a bipolc 
will tend to favour the hemispheric pattern. Importantly, however, the bipolc 
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Figure 5. Shape of the field lines around the PIL of a bipole in a background field. Colours 
show the skew angle along the PIL for each possible sign combination of bipole twist and 
background field polarity. The colour table runs from red (dextral skew) through green (weak 
skew) to blue (sinistral skew). In panels a and d (left column) the bipole has positive twist /3, 
in panels b and e (centre column) the bipole is untwisted (/3 = 0), and in panels c and f (right 
column) the bipole has negative twist. In panels a, b, and c (top row) the following polarity 
of the bipole matches that of the background field, and in panels d, e, and f (bottom row) the 
leading polarity matches the background. The shading shows polarity of radial magnetic field 
on the photosphere (white positive, grey negative), and thin lines show contours of the same 
field. Thick lines are selected coronal magnetic field lines crossing the PIL. 

shape can produce exceptions to the hemispheric pattern for filaments occuring 
at either the north or south end of a bipole, depending on the relative polarity 
of the background field. 

3.4. Multiple Bipole Interactions 

The shape of a single region has a strong influence on the skew of the field 
above its external PIL, as described in Section 3.3. However, in a minority 
of cases a non-uniform background field modifies how the bipole connects to 
other nearby bipoles or regions of strong flux. This is mechanism 5 in Figure 3, 
where the example shows a bipole which has emerged to the north of a strong 
negative region. Rather than connecting directly across the PIL, the field lines 
from the positive leading polarity have been diverted southward to connect with 
the stronger region. Thus a part of the PIL which would normally be expected 
from Figure 5 to have weak skew actually has strong dextral skew - in fact this 
was the location of an observed dextral filament. Of course, this mechanism 
may also produce sinistral skew, as would have been the case had the strong 
background region been located to the north of the new bipole. 

In addition to the shape of a bipole being modified by the pre-existing field, 
the skew of the pre-existing field may also be modified by a new bipole emergence, 
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through the dynamic effect of sweeping (mechanism 6, Figure 3). The technique 
of sweeping before the insertion of a new bipole was described in Papers 1 and 
2, and results in a more realistic end configuration for newly-emerged regions. 
In the example shown here in Figure 3, a new region emerges to the left of 
a pre-existing PIL with slight dcxtral skew. The pre-existing positive field has 
been swept out of the insertion region, resulting in (1) convergence of footpoints 
towards the pre-existing PIL, and (2) lengthening of field lines across this PIL 
because the positive footpoints are displaced northward. The result is a highly 
skewed dcxtral flux rope structure along the PIL. Again, this mechanism could 
also produce the opposite sign of skew, and is also capable of introducing skew 
at a PIL where there was previously none. Although our sweeping mechanism 
is artificial, we might expect a similar effect on pre-existing coronal field even 
from a gradual flux emergence. 

It was described in Section 3.1 how differential rotation is the primary surface 
transport effect that develops skewed coronal fields from scratch, and that surface 
diffusion acts primarily to enhance existing skew. In fact, in a small number of 
cases, diffusion does cause skew to develop from scratch, and this is mechanism 
7 of Figure 3. The necessary ingredients are a PIL with initially complex shape. 
Diffusion inevitably smooths out the PIL, but the coronal field is not necessarily 
smoothed to the same extent. Thus the final PIL shape may end up sheared 
with respect to the overlying coronal field lines. Wang and Muglach (2007) show 
an observed sequence in their Figures 2 and 3 where this mechanism appears to 
occur. 

3.5. Skew from the Initial Potential Field 

Finally, mechanism 8 in Figure 3 relates to the initial condition of the simulation, 
which is a potential field extrapolation. Although this field has the minimum en- 
ergy for the boundary conditions, and so the majority of field lines are unsheared 
and pass over PILs at 90°, this is not necessarily always the case. The example 
in Figure 3 shows an example of how skew may arise in such a potential field. 
Although the lower field lines are not skewed, higher coronal field lines must link 
regions of strong flux, so in this case are skewed with respect to the PIL because 
of the relative locations of the two strong-flux regions. On the real Sun — with 
no initial potential field — it is unclear whether these overlying field lines would 
follow the same path. 

4. Relative Importance of Different Mechanisms 

To assess the role of each of the eight mechanisms in producing the hemispheric 
pattern, we present here the results of a detailed analysis of the 109 filaments 
which were compared to observations in Paper 2. Table 1 shows the number 
of filaments in relation to each mechanism. Note that the mechanisms are not 
mutually exclusive; there are 56 filaments with two or more mechanisms associ- 
ated, and 12 with four or more. The table also breaks various statistics down by 
mechanism, and these are described in this section. The bottom row shows the 
same statistics for the full sample of 109 filaments. 
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Table 1. Classification of filaments by skew mechanism (as defined in Figures 2 and 3). Sim- 
ulation results presented here are for run 4 (bipole twist = 0.4, with observed hemispheric 
sign of helicity), except for the fifth column which compares run 4 with run 1 (|/3| = 0.2, 
opposite sign of helicity). 
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1.85 
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109 
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32.1 


1.58 



The most common mechanism found is mechanism 4 (bipole shape), followed 
by mechanisms 3 (internal filaments) and 1 (differential rotation at a north- 
south PIL). By contrast, differential rotation at an east-west PIL (mechanism 2) 
is relatively rare, as is the creation of skew by surface diffusion of a complex PIL 
(mechanism 7). Note that mechanism 2 may be slightly more common on the 
real Sun since our sample of 109 filaments is biased against those on the polar 
crown, due to projection effects hindering chirality observations. 

The third column of Table 1 shows the percentage of simulated filaments 
following the hemispheric pattern. This is for the simulation run with |/3| = 0.4 
and the observed hemispheric sign of bipole helicity (run 4 in Paper 2). This was 
the run which best agreed with the observed filament chiralities. As indicated in 
Section 3, mechanisms 1 and 3 always produce skew following the hemispheric 
pattern, and mechanism 2 produces exceptions. The percentage for mechanism 
2 is above zero because it was not the dominant mechanism in all cases where 
there were multiple factors at play. As expected, mechanisms 4 to 8 produce some 
filaments following the hemispheric pattern and some not. However, we see that 
the percentages of filaments following the pattern are in the range 70%-88%, 
which explains the overall result. The lowest percentage for these mechanisms 
is mechanism 4 (bipole shape), and it was shown in Section 3.3 how this can 
naturally be responsible for exceptions to the hemispheric pattern. The question 
arises as to how mechanisms 5, 6, and 7 preferentially give the correct pattern, 
even though they can produce both types of skew. The answer is that they are 
usually not the only mechanism acting. They tend to enhance existing skew 
rather than generate completely new skew from scratch. Thus mechanisms 1, 3, 
and 4 (differential rotation and bipole shape) arc more fundamental. 

Column 4 shows how the simulated filaments involving each mechanism match 
with the observed chirality at those locations. The agreement is generally good 



arymechanisms_rev.tex; 2/10/2008; 19:17; p. 10 



Modelling the Global Solar Corona III 



(above 85%). Filaments with mechanism 2 (east- west PIL) are least likely to 
be correct, supporting the idea that this mechanism is not responsible for the 
chirality of filaments on the Sun. In the simulation this mechanism is associated 
primarily with high-latitude filaments forming from the initial potential field, 
whereas in real-life we expect that the correct hemispheric sign of skew will be 
transported polewards from lower latitudes over longer periods of time. 

The fifth column of Table 1 shows the relative importance of emerging active- 
region hclicity for the filament hemispheric pattern. It compares a run with 
observed hemispheric sign of (3 (run 4 in Paper 2) against a run with the opposite 
sign of (3 (run 1). A higher percentage means that there are more filaments where 
changing the bipole twist j3 has an effect on the resulting chirality. From the 
column we see that (3 has the most influence on internal filaments (mechanism 
3), and on filaments formed by mechanisms 4, 5, and 6. These mechanisms 
are all associated with the emergence of new bipolar regions. However, even 
for mechanism 3 the percentage is only 51.4%, because the effect of emerging 
twist is over-ridden after 1 to 2 Carrington rotations by differential rotation (as 
described in Section 3.2). This shows the importance of carrying out continuous 
long-term simulations rather than individual extrapolations. 

For mechanisms 2,7, and 8 the emerging twist has less influence, as these are 
largely associated with the initial potential field. Overall only 32% of filaments 
change their skew type when the hemispheric sign of emerging bipole twist is 
reversed, and even if filaments related to the initial field are discounted, the 
percentage will still be below fifty. Thus the emergence of twisted magnetic 
fields is not the most important factor in producing the hemispheric pattern, 
especially for the majority of filaments which form between multiple bipolar 
regions. This explains why our simulation with the wrong sign of (3 (run 1 in 
Paper 2) still produced the correct chirality for over 60% of filaments. 

Each filament in the sample was associated with a subset of the 119 emerging 
bipolcs in the simulation. These were selected as being those which influenced 
the skew above the PIL in some way, rather than simply being located nearby. 
The final column of Table 1 gives the mean number of bipoles associated with 
each filament, again broken down by mechanism. Not surprisingly, mechanisms 5 
(relative positions of flux) and 6 (sweeping) are associated with the most bipolcs, 
while mechanisms 2 and 3 have the fewest. In the case of mechanism 2 (east- west 
PIL) this is because such locations are at high latitudes and created from the 
initial potential field only. The overall mean is 1.58 regions, demonstrating that 
most filaments are associated with more than one active region. This supports 
the result of Mackay et al. (2008). 

Finally, consider how Table 1 would differ if (3 = 0, i.e. if the bipolcs emerge 
untwisted (run 2 in Paper 2) . The same mechanisms generally occur in that run, 
except that mechanism 3 (internal PIL) generates skew more slowly because of 
the lack of initial twist. Indeed, this explains 7 out of the 10 filaments found 
to be weakly skewed in run 2 (sec Table 3 in Paper 2). As a consequence, the 
overall percentages of filaments following the hemispheric pattern and matching 
the observed chirality would decrease, but only slightly. Figures 5(b) and (e) 
show that, even for an untwisted bipole, there is still skew at either end of the 
external PIL, so that mechanism 4 (bipole shape) and other mechanisms still 
occur. 
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5. Conclusion 

In this paper we have presented the results of a detailed analysis of how skew 
is formed at the filament locations in our 3D simulation of the global solar 
coronal evolution over 6 months (see Paper 2 for details of the simulation). The 
simulation produced the correct chirality for up to 96% of the 109 filaments 
compared with observations, an unprecedented result. We have identified eight 
key mechanisms that produce skewed coronal fields above PILs, and these are 
summarised in Figures 2 and 3. The relative importance of these mechanisms 
has been assessed and we are now in a position to explain why the simulation 
produces the correct hemispheric pattern. With the exception of mechanism 
8 (connected to our initial potential field), we suggest that all of these mech- 
anisms occur in some form on the real Sun, so that our explanation for the 
simulated hemispheric pattern should offer some insight into the origin of the 
real hemispheric pattern. 

There is not one single factor responsible for the pattern, but rather a number 
of factors depending on the location of each individual filament and the nature 
of the PIL. The pattern arises as a statistical result because of the relative 
frequency of these mechanisms. The most important single influence is the bipo- 
lar shape of active regions, which is a well-observed feature, except for some 
activity complexes. The location where a filament forms around or inside such a 
bipole has a key influence on its skew type. We have shown that on average the 
hemispheric pattern is upheld for filaments forming either on the internal PIL of 
a bipole, or around its external edge. The second key mechanism is differential 
rotation, but acting on a north-south PIL rather than an east-west PIL, as 
the latter would tend to produce exceptions to the hemispheric pattern. For 
a north-south PIL, differential rotation leads naturally to dextral skew in the 
northern hemisphere and sinistral in the southern hemisphere. Such north-south 
PILs occur frequently both within individual bipoles and between neighbouring 
active regions, the locations where many filaments are observed to form. 

Once an initial skew is produced by one of the fundamental mechanisms 
described above, it is often enhanced by the action of surface (supergranular) 
diffusion, which converges the footpoints of coronal field lines towards the PIL, 
increasing the skew angle. Depending on individual circumstances, the skew may 
also be either enhanced or weakened by one of the remaining mechanisms in our 
list, such as the distribution of strong flux regions in the background field, or 
nearby bipole insertions. In some cases the lift-off of a twisted magnetic flux rope 
(as in Mackay and van Ballcgooijcn, 2006) was found to remove skew, though 
typically only from higher field lines. It was found that the original source of skew 
at each filament location could develop many days before the filament's date of 
observation (over 100 days for many filaments in this simulation). Over this time 
a number of emerging active regions could affect the skew of the field through 
the various mechanisms described. This highlights the importance of simulating 
the long-term coronal evolution, rather than doing single-time extrapolations. 

A fundamental problem which must be addressed by any theory hoping to 
explain the hemispheric pattern is how exceptions to the pattern may be formed. 
The results of our study demonstrate that in our model the pattern holds in 
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a statistical sense, and that exceptions are a natural feature. Out of the 109 
filaments in our sample, there were 16 with observed chirality opposite to the 
hemispheric pattern. The majority of these formed on the outside edge of a 
bipolar region, and six are simply explained by the shape of the bipolar field, as 
explained in Section 3.3. Three others involve both the shape of the bipole and 
its relation to the surrounding field (mechanisms 4, 5, and 6 in Figure 3), giving 
examples of how such interactions can produce exceptions to the hemispheric 
pattern in some circumstances. Of the remaining filaments, four involve flux from 
the initial potential field, where our simulation is not expected to be realistic, 
and two are internal to individual bipoles. We assume here that on the real 
Sun these two active regions must have had the minority sign of helicity for 
their hemisphere; in the simulation run with opposite sign of bipole twist we 
reproduce the observed chiralitics for these filaments. 

In summary, we have demonstrated how the hemispheric pattern of filament 
chirality on the Sun may be explained by surface flux transport and the emer- 
gence of twisted active regions. The key elements are (1) differential rotation 
at north-south PILs and (2) the shape of bipolar active regions, along with an 
evolution over periods of many days. However, there were many complex exam- 
ples where interactions between a number of active regions became important. 
These interactions involved multiple bipoles of different sizes, relative locations, 
and dates of emergence, including but not limited to the particular two-bipolc 
configurations suggested by Martens and Zwaan (2001). Our simulations support 
their view (and that of Zirker et al., 1997) that the skew above PILs on the Sun 
develops over a period of months, through interaction and merging of multiple 
active regions. There are observed cases where strong skew appears to develop 
within a day or two (e.g. Gaizauskas et al, 1997), and this sometimes occurs in 
the simulation following new bipole emergence. However, by comparing runs with 
opposite signs of emerging bipole helicity, we find that the sign of (3 affects the 
chirality of only about one third of filaments, conflicting with simple emerging 
flux rope explanations for the hemispheric pattern (Rust and Kumar, 1994). 

Finally, we note that our 6-month simulation does not produce the correct 
chirality for high-latitude filaments, including those on the polar crown. In the 
present simulation, differential rotation of the initially potential field produces 
the opposite skew across these primarily east-west PILs. This opposite sign of 
skew at high latitudes is also reflected in a sign-reversal of current helicity at 
these latitudes, which has been found to persist for over two years in longer 
simulations (Ycates, Mackay, and van Ballegooijcn, 2008b). Interestingly, obser- 
vations of coronal X-ray arcades by McAllister et al. (1998) found their skew 
to agree with our present simulations, and to be opposite to their underlying 
filament chiralitics. Future simulations will need to consider whether poleward 
helicity transport over the full solar cycle is enough to counteract the effect of 
differential rotation on the chirality of high-latitude filaments. 
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